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Measurements of heavy vector boson pair production 
(WW, WZ, ZZ) are of great importance because they 
test the electroweak sector of the standard model (SM). 
Diboson production provides a sensitive probe of new 
physics, including anomalous trilinear gauge couplings 
new particles such as the Higgs boson [2J], and large 
extra dimensions 0]. Important tests of electroweak 
physics have been made recently at the Fermilab Teva- 
tron with the first observation of WW production in 
hadron collisions 0, @| and the first observation of WZ 
production The production of Z pairs has been ob- 
served in e + e~ collisions at LEP Q, but not in hadron 
collisions. As a window to new physics, ZZ production 
is particularly interesting because of the absence of ZZ^f 
and ZZZ couplings in the SM, and because of the very 
low backgrounds in the four charged-lepton channel. 

The most sensitive previous search for ZZ production 
in hadron collisions was reported by the DO Collabora- 
tion using data corresponding to 1 fb _1 of integrated lu- 
minosity. That search used only the four charged-lepton 
channel and set a limit on the cross section of a(ZZ) < 
4.4 pb at 95% C.L. 8]. The next-to-leading order (NLO) 
ZZ cross section for pp collisions at y/s = 1.96 TeV is 
1.4 ± 0.1 pb in the zero-width Z approximation 

In this Letter, we report a signal for ZZ production in 
hadron collisions which has a probability of 5.1 x 10~ 6 
to be due to the expected background. The production 
is observed in pp collisions at i/s = 1.96 TeV using a 
data sample corresponding to 1.9 fb _1 of integrated lu- 
minosity collected by the CDF II detector. We consider 
both ZZ -> £££'£' and ZZ -► llvv decays, where t and £' 
are electrons or muons directly from Z decay or from the 
leptonic decay of r's when one or both Z bosons decay 
to t leptons. 

The CDF II detector [l(| geometry is described us- 
ing the azimuthal angle <f> and the pseudo-rapidity 77 = 
— ln[tan(0/2)], where is the polar angle with respect 
to the proton beam axis (positive z-axis). The pseudo- 
rapidity of a particle assumed to originate from the center 
of the detector is referred to as 

The trajectories of charged particles (tracks) are recon- 
structed using silicon micro-strip detectors fill . [l2T | and 
a 96-layer open-cell drift chamber (COT) [l3| inside a 
1.4 T solenoid. A particle with < 1 traverses all 96 
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COT radial layers. The number of traversed layers de- 
creases to zero as \r)&\ approaches 2. The silicon system 
provides coverage with 6 (7) layers with radii between 2.4 
cm and 28 cm for \r] d \ < 1.0 (1.0 < |r?d| < 2.0). Outside 
of the solenoid are electromagnetic (EM) and hadronic 
(HAD) sampling calorimeters, segmented in a projective 
tower geometry, and constructed of layers of lead and 
iron absorber, respectively, and scintillator. The EM sec- 
tion is the first 19-21 radiation lengths, corresponding to 
one hadronic interaction length (A), and contains elec- 
tromagnetic showers. The HAD section of the calorime- 
ter extends to 4.5-7 A and contains the majority of a 
hadronic shower. The calorimeters are divided into cen- 
tral (\r)d\ < 1.1) and forward (1.1 < |?y d | < 3.6) regions. 
Outside of the central calorimeters are muon detectors 
consisting of scintillators and drift chambers. 

The branching fraction of the ZZ state to four e or fi 
leptons, including those from leptonic t decays, is only 
0.51%. When coupled with the small SM cross-section, 
only a small number (^14) of ZZ — > £££'£' are expected 
to be produced in 1.9 fb _1 at the Tevatron. The finite 
CDF II acceptance further reduces the expected number 
of observed ZZ ->■ £££'£' events. 

To maximize the acceptance, we construct lepton can- 
didates out of all reconstructed tracks and energy clusters 
in the EM section of the calorimeter. This is done with 
the same lepton identification criteria used in a previous 
CDF measurement of WZ production 0. The lepton 
candidates are divided into seven exclusive categories: 
three for muons, three for electrons, and one for "track- 
only leptons" which are tracks extrapolating to detec- 
tor regions that are inactive for energy measurement be- 
cause they are either not covered, or are only partially 
covered, by calorimeter components. The three muon 
categories include one that uses the muon detectors and 
two that use different types of minimum ionizing tracks, 
central (\r] d \ < 1.1) and forward (1.2 < |?7 d | < 2.0). The 
three electron categories are central ( 1 77,^ | < 1.1), for- 
ward (1.2 < |77d| < 2.0) with a matched silicon-based 
track, and forward (1.2 < |??dj< 2.8) without a matched 
silicon-based track (see Ref. [6J for more detail). 

The transverse energy Et of an EM cluster or 
calorimeter tower is Esm6, where E is the associated 
energy. Similarly, px is the component of track mo- 
mentum transverse to the beam line. To suppress jets 
being misidentified as leptons, all lepton candidates are 
required to pass a calorimeter-based isolation require- 
ment such that the sum of the Et for the calorime- 
ter towers not associated with the lepton in a cone of 
AR = y/ (Arf) 2 + (Ac/)) 2 < 0.4 around the lepton direc- 
tion is less than 10% of the Et for electrons or pt for 
muons and track lepton candidates. 

To detect the presence of neutrinos in ZZ — > ££vv de- 
cay, we use missing transverse energy $J T = — ^T,j, 
where fiT.i is the transverse component of the unit vector 
pointing from the interaction point to calorimeter tower i. 
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The ]firp calculation is corrected for muons and track-only 
lepton candidates, which do not deposit all of their en- 
ergy in the calorimeter. We analogously define the scalar 
sum T,E T = JV E T ^ applying the same corrections. 

The events we consider must pass at least one of four 
trigger selection criteria. The central electron trigger re- 
quires an EM energy cluster with Et > 18 GeV matched 
to a track with pt > 8 GeV/c. For the ZZ — ► eevv chan- 
nel only, a trigger for forward electrons requires an EM 
energy cluster with Et > 20 GeV and an uncorrected, 
calorimeter-based measurement of $ T > 15 GeV. Two 
muon triggers are based on stubs from the corresponding 
muon detectors matched to a track with px > 18 GeV/c. 
Trigger efficiencies are measured in leptonic W and Z 
data samples [14 1. 

The ZZ — > £££'£' candidates are selected from events 
with exactly four charged-lepton candidates using re- 
quirements that were optimized with Monte Carlo sim- 
ulation without reference to the data. At least one lep- 
ton is required to satisfy the trigger criteria and have 
E T > 20 GeV {p T > 20 GeV/c) for electrons (muons). 
We loosen this requirement to 10 GeV (GeV/c) for the 
other leptons to increase the ZZ kinematic acceptance. 
We require at least two same-flavor, opposite-sign lepton 
pairs in the event. Trackless electrons are considered to 
have either charge, and track-only leptons either flavor. 
One pair must have invariant mass Mg+i- in the range 
[76, 106] GeV/c 2 , while the requirement for the other 
pair is extended to [40, 140] GeV/c 2 to increase the ac- 
ceptance for off-shell Z decays. 

The acceptance for the ZZ process is determined us- 
ing PYTHIA 15j Monte Carlo calculations followed by a 
GEANT-based simulation [l6| of the CDF II detector. An 
efficiency correction, of up to 10% per lepton, is applied 
to the simulation based on measurements of the lepton re- 
construction and identification efficiencies using observed 
Z -> £ + £~ events. 

The dominant backgrounds to the ZZ — > £££'£' selec- 
tion are the Drell-Yan Z/j* process (DY) with two jets 
misidentified as leptons (Z+jets) and DY with an ad- 
ditional photon and a jet, both misidentified as leptons 
(Z7+jets). The Z+jets and -Z^+jets background contri- 
butions are estimated from data by extrapolating from 
a sample of events that contain three identified leptons 
and a jet ji containing a track or EM energy cluster sim- 
ilar to those required in the lepton identification. The 
contribution of each event to the total yield is scaled by 
the probability that the ji is identified as a lepton. This 
probability p(ji) is determined from multijet events col- 
lected with jet-based triggers and is a function of the ji 
Pt and type of lepton. A correction to p{ji) is applied 
for the small real lepton contribution using Monte Carlo 
simulation of single W and Z boson processes. In this 
background sample, one of the three identified leptons is 
likely to be either a jet or a photon misidentified as a 
lepton. An event with two leptons and two j\ jets enters 



TABLE I: Expected and observed number of ZZ -> Llt'f 
candidate events. The first uncertainty is statistical and the 
second one is systematic. 





Candidates without a 


Candidates with a 


Category 


trackless electron 


trackless electron 


ZZ 


1.990 ± 0.013 ± 0.210 


0.278 ± 0.005 ± 0.029 


Z+jets 


0.014tom° ± 0.003 


0.082±o;°60 ± °- 016 


Total 


2.004±°°^ ± 0.210 


0.360tg;gio ± 0.033 


Observed 


2 


1 



the three lepton plus ji sample if either of the ji jets is 
misidentified as a lepton, but will enter the four identified 
lepton sample only if both ji jets are misidentified as lep- 
tons. Therefore, the contribution from this category of 
events is double counted. A correction for this is made by 
subtracting the yield of two leptons plus two ji jets scaled 
by p(Ji,i) Xp(jl,2)- As the three lepton plus ji sample has 
significant contributions from the ZZ — > £££'£' signal it- 
self when one of the leptons is not fully identified but is 
counted as a ji, an anti-isolation requirement of >20% 
(see previous definition of calorimeter-based isolation) is 
applied to the ji selection. 

As a cross-check, we estimate the Zj+jet background 
contribution in an alternative way using the yield of three 
lepton plus ji events in simulated Zj+jet data scaled by 
p(ji). We correct this with the ratio of data to simulation 
for an analogous calculation of two identified leptons and 
one ji scaled by p(ji) in Z+jets events. This estimate of 
the Zj+jet background is in good agreement with the 
nominal estimate based solely on the data. 

We separate the ZZ — » £££'£' candidates into two 
exclusive categories based on whether or not they con- 
tain at least one forward electron without a track. This 
is done because the background from Z^+]ets is much 
larger in candidates with a forward trackless electron. 
The expected signal and background yields assuming 
cr(ZZ) = 1.4 ± 0.1 pb and the observed yields are shown 
in Table [J 

The statistical significance of the ZZ — > £££'£' yield is 
determined using a maximum likelihood fit with two bins, 
one for each of the ZZ — > £££'£' categories. We define 
Aln£ as the logarithm of the likelihood ratio between 
this fit and the no signal hypothesis. In 10 7 background- 
only Monte Carlo experiments, only 109 have larger 
A In £ than that observed in data. This corresponds to a 
background-only probability (p- value) ofl.lxlO -5 and a 
signal significance equivalent to 4.2 standard deviations. 

The ZZ — * ££vv candidates are selected from events 
with exactly two lepton candidates excluding events with 
forward electrons without a track which are contami- 
nated by large IV7 backgrounds. At least one lepton 
is required to satisfy the trigger and have Et > 20 GeV 
(pt > 20 GeV/c) for electrons (muons). This require- 
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merit is loosened to 10 GeV (GeV/c) for the other lep- 
ton. We apply a track-based isolation selection in which 
the sum of the px of the tracks not associated with the 
lepton within a cone of AR < 0.4 around the lepton is 
required to be less than 10% of the momentum of the 
track associated with the lepton. 

Aside from ZZ production, other SM processes that 
can lead to two high-p^ leptons include events from DY, a 
W decay with photon (W-f) or jet (W+jets) misidentified 
as a lepton; and it, WW, and WZ production. The it 
contribution is suppressed by requiring fewer than two 
reconstructed jets with Et > 15 GeV and \r]d\ < 2.5 in 
the event. The DY background is suppressed by requiring 
sufficiently large $! T in the event to remove contributions 
from mismeasured leptons and/or jets. This is achieved 
by requiring $ T rel > 25 GeV, where 



$ T sin 



,(^,jet) 



6r,(^Jet) > f 
5 T ,(£,jet) < f 



(1) 



and A(^>* t (gj e t) 1S the angle between the $ T direction 
and the nearest lepton or jet. To suppress events with 
Ifjrp from mismeasured unclustered energy, we require sig- 
nificant Jfjrp such that $ T l^f{T 1 E T ) > 2.5 GeV 1 / 2 . We 
require the lepton pair to be consistent with the same- 
flavor, opposite-sign property of Z decay and have dilep- 
ton invariant mass M^+i- > 16 GeV/c 2 to suppress QCD 
backgrounds. 

The acceptances for the WW, WZ, ZZ, W~f, and it 
processes are determined using the same detector simu- 
lation as described for the ZZ — » III' I' channel. Events 



are simulated with the MC@NLO program for WW [17j |. 
PYTHIA for WZ, ZZ, and it [HI, and the generator de- 
scribed in [li| for W'y. An additional correction is ap- 
plied to the W-f background estimate based on a mea- 
surement of the photon conversion veto efficiency in data. 
The background from FK+jets is estimated from the yield 
of one identified lepton plus one ji scaled by p(ji). As the 
sample size is sufficiently large, a loose calorimeter-based 
isolation cut (< 30%) is applied to the ji samples to re- 
duce the magnitude of the extrapolation from the ji to 
the fully-identified lepton. 

We observe 276 events in the selected region which 
is expected to contain 256 ± 21 events of which only 
14 ± 2 are from the ZZ — > ilvv process in the SM. Ap- 
proximately half of the yield is due to the WW process. 
However, ZZ — + llvv and WW have different kinematic 
properties which are exploited to statistically separate 
the contribution of these two processes to the dataset. 
We calculate an event-by-event probability density for 
the observed lepton momenta and $ T using leading order 
calculations of the differential decay rate for the processes 
The event probability density is 



P(x obs ) = 



dcr L o{y) 
dy 



S3 10 
c 

CD 



•Data 0WW f L dt = 1.9 fb 1 

□ZZ ffltt,Wy,W+jets 
□WZ 0DY 

HSyst. Uncertainty 




-5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 

log(l-LR) 

FIG. 1: Distribution of the discriminating variable log 10 (l 
LR) for the ZZ — > llvv search. 



where the elements of y (x fcs) are the true (observed) 



values of the lepton momenta and $ T . 



d<T L 



dy 



is the par- 



e(y)G(x obs ,y)dy, (2) 



ton level cross section differential in those observables, 
e(y) is the detector acceptance and efficiency function, 
and G(x bs,y) is the transfer function representing the 
detector resolution. The constant (a) normalizes the to- 
tal event probability to unity. The missing information 
due to the fact that we have two neutrinos in the final 
state is integrated over in this calculation. We then form 
a likelihood ratio discriminant LR which is the signal 
probability divided by the sum of signal and background 
probabilities LR = Pzz/(Pzz + P\vw)- The distribution 
of log 10 (l — LR) for the data compared to the summed 
signal and background expectation is shown in Figure [T] 

For both III' I' and llvv channels, the systematic un- 
certainties associated with the Monte Carlo simulation 
affect the WW, WZ, ZZ, W-f, and it expectations sim- 
ilarly. The uncertainties from the lepton selection and 
trigger efficiency measurements are propagated through 
the analysis, giving uncertainties from 1.0% to 1.4% and 
2.1% to 6.1% for the respective efficiencies of the differ- 
ent signal and background processes. The detector ac- 
ceptance variation due to PDF uncertainties is assessed 
to be in the range 1.9% - 4.1% using the 20 pairs of PDF 
sets described in We assign a 6% luminosity uncer- 
tainty to signal and background estimates obtained from 
simulation [20j. The cross-section uncertainties are 10% 
for WW 0, WZ @, and W-f [Hj], and 15% for it [HHH. 

The systematic uncertainty on the VF+jets background 
to ZZ — > llvv candidates and the Z+jets and Z7+jets 
background to the ZZ III' I' is estimated to be 20% 
from differences in the observed probability that a jet 
is identified as a lepton for jets collected using different 
jet Et trigger thresholds. These variations correspond 
to changing the parton composition of the jets and the 
relative amount of contamination from real leptons. 

For the ZZ — > llvv, the systematic uncertainty on 
the DY background yield due to the $ T resolution mod- 
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eling is estimated to be 20% from comparisons of the 
data and Monte Carlo simulation in a sample of dilepton 
events. For the Wj background contribution, there is an 
additional uncertainty of 20% from the detector material 
description and conversion veto efficiency. 

We define four independent control samples based on 
the ZZ — > llvv selection where one of the cuts is re- 
moved or inverted to test our modeling of background- 
dominated data. Removing the $ T requirement produces 
a DY-dominated sample which tests luminosity account- 
ing, lepton reconstruction efficiency, and non-^ T related 
trigger efficiencies that apply to both the III 1 1' and llvv 
final states; we observe (expect) 160980 (160000± 18000) 
events, where the uncertainty combines statistical and 
systematic contributions. Inverting the charge sign re- 
quirement to select same-sign dilepton events tests our 
modeling of photons and jets misidentified as leptons; we 
observe (expect) 161 (138 ± 19) events. Inverting the 
$ T / y^'EEr) requirement but selecting $ T rel > 25 GeV 
events tests our modeling of the effect of unclustered en- 
ergy on the $ T ; we observe (expect) 55 (59 ± 9) events. 
Finally, inverting the lfi T rel > 25 GeV requirement but 
selecting $ T > 25 GeV events tests our modeling of mis- 
measured leptons or jets at high $ T ; we observe (expect) 
151 (178 ± 30) events. The observed yields are in good 
agreement with the expectations in each selection set. 

We combine the ZZ -> llvv with the ZZ -> III' I' re- 
sults by extending the likelihood fit previously described 
to include the log 10 (l — LR) distribution. The p- value 
for the ZZ — > llvv alone is 0.12 and the combined p- 
value is 5.1 x 10~ 6 corresponding to a significance equiv- 
alent to 4.4 standard deviations. We determine a ZZ 
cross-section for this signal by fitting the data for the 
fraction of the expected SM yield in the full acceptance 
and scaling the zero- width Z boson approximation cross- 
section by that fraction. The measured cross section is 
<j{pp — ^ ZZ) = 1.4^ 6 pb, consistent with the SM ex- 
pectation. This is the first measurement of the ZZ cross 
section in hadron collisions. 
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